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Abstract 

We present a new mechanism for generating the baryon asymmetry of the universe 
directly in the decay of a singlet scalar field S r with a weak scale mass and a high 
dimensional baryon number violating coupling. Unlike most currently popular models, 
this mechanism, which becomes effective after the electroweak phase transition, does 
not rely on the sphalerons for inducing a nonzero baryon number. CP asymmetry in 
S r decay arises through loop diagrams involving the exchange of gauge bosons, 
and is suppressed by light quark masses, leading naturally to a value of tjb ~ 1CP 10 . 
We show that the simplest realization of this mechanism, which uses a six quark AB = 
2 operator, predicts colored scalars accessible to the LHC, and neutron-antineutron 
oscillation within reach of the next generation experiments. 



1 Introduction 



Recent developments in particle physics have had profound impact on cosmology. One of 
the most far-reaching consequences has been the possibility that new interactions beyond 
the standard model can explain the origin of matter-antimatter asymmetry of the universe 
as a dynamical phenomenon. There are currently several attractive scenarios which achieve 
this, the two most widely discussed ones being (i) baryogenesis via leptogenesis [I], which is 
connected to the seesaw mechanism and neutrino masses, and (ii) weak scale baryogenesis 
which involves supersymmetric or multi-Higgs extensions of the standard model. Both 
these proposals depend crucially on the properties of the elctroweak sphaleron [3| which 
serves as the source of B violation. Since the nature of new physics beyond the standard 
model remains unknown presently, it is important to explore alternative mechanisms that 
can explain the matter-antimatter asymmetry while yielding testable consequences. In this 
letter we suggest and explore one such alternative. 

The salient feature of our proposal is that baryogenesis occurs via the direct decay of a 
scalar boson S r having a weak scale mass and a high dimensional baryon violating coupling. 
S r is the real part of a baryon number carrying complex scalar S, which acquires a vacuum 
expectation value (vev). The decays S r — > 6q and S r — > 6q will then be allowed, providing 
the source for B asymmetry. These decays occur when the temperature of the universe 
is T ~ 0.1 — 100 GeV. By this time the electroweak sphalerons have gone out of thermal 
equilibrium, and thus play no role in the B asymmetry generation. We call this mechanism 
"post-sphaleron baryogenesis" . The three Sakharov conditions for successful baryogenesis jlj 
are satisfied rather easily in our scheme. The high dimensionality of the B violating coupling 
of S r to the quark fields allows the AB ^ decays to go out of equilibrium at weak scale 
temperatures. CP violation occurs in the decay via loop diagrams involving the exchange 
of the standard model W ± gauge bosons. This amplitude has sufficient light quark mass 
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suppression to explain naturally the observed (small) value of the baryon to photon ratio 
i]b ~ 10~ 10 . The simplest realization of our mechanism involves interactions that violate B 
by two units and therefore gives rise to neutron-antineutron oscillations. We find that the 
successful implementation of our mechanism sets an upper limit on the transition time for 
N <r-> N oscillation bringing it to within the realm of observability This connection provides 
a strong motivation for improved searches for iV <-> N oscillation [3] . 

The connection with N <-> N oscillation can be understood as follows. Let us consider an 
interaction of the form SO^b, where S is a standard model singlet complex scalar field (with 
S r denoting its real part) and Oab is the baryon number violating operator in question. This 
interaction will lead to baryon number violation if (S) ^ 0. Suppose the mass dimension of 
the operator SO^b is M~ n with n positive. The higher the value of n for an operator O^b, 
the lower the mass scale allowed by the existing limits on baryon violation. Since the rate 
of these AB ^ interactions in early universe goes like M~ 2n , the higher the value of n, 
the easier it is to satisfy the out-of-equilibrium condition at a lower temperature (multi GeV 
range). Clearly the operator leading to B — L conserving proton decay mode cannot be useful 
for us, since present experimental limits on proton lifetime imply that this operator should 
go out of equilibrium at temperatures of order 10 14 — 10 15 GeV. On the other hand, for a 
process like N <->• N oscillation jHl 13 E] , present experimental lower limits on the oscillation 
time r N _fj allow the mass M appearing in the operator (f gh)u c d c d c u c d c d c / M 5 to be 

in the multi-TeV range (for the first family Yukawa couplings / ~ g ~ h ~ 10~ 3 ). The 
out-of-equilibrium temperature for the processes S r — > 6q and S r — > 6q is then allowed to 
be below the sphaleron decoupling temperature of about 100 GeV. We will illustrate how 
post-sphaleron baryogenesis works using the AB = 2 process, although the mechanism 
applies more generally. (The AB = 2 operator involving left-handed quark doublet fields, 
QQQQQQH* H* /M 7 , has additional Higgs fields and thus a higher dimensionality.) 
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The high dimensional AB = 2 couplings of S are obtained by integrating out colored 
scalar fields. These colored scalars cannot be much heavier than about a TeV, or else the 
induced t]b, consistent with nucleosynthesis limits, will be too small. The prospects for 
discovering such baryon number carrying colored scalars at the Large Hadron Collider are 
quite promising. 

An attempt to generate baryon asymmetry at a temperature of order MeV via the decay 
of a heavy (~ 50 TeV) gravitino within supergravity was proposed in Ref. jllj . Such a large 
gravitino mass would however require fine-tuning to solve the hierarchy problem. Another 
scenario ^2] invokes the decay of the inflaton into squarks, with their subsequent decay 
producing baryon asymmetry. This mechanism requires that the reheating temperature be 
less than a GeV in order for the scattering and inverse decays not to wash out the asymmetry. 
The model presented here differs from these earlier attempts in two crucial ways: (i) There is 
a strong link between baryon asymmetry and N <-> N oscillation, and (ii) the mechanism of 
inducing CP asymmetry via the standard model loops which leads naturally to a value 
of r\B ~ 10~ 10 due to light quark mass suppression is entirely new. 

2 Light diquarks and observable N <-> N oscillation 

To illustrate our mechanism for post-sphaleron baryogenesis, we consider a generic TeV 

scale model that gives rise to the higher dimensional operator for N <-> N oscillation. It 

consists of the following color sextet, SU(2)l singlet scalar bosons (X, Y, Z) with hypercharge 

— §,+§>+§ respectively that couple to the right-handed quarks. 1 In addition, there is a 

complex scalar field S which is a singlet of the standard mode with mass in the 100 GeV 

range. With this field content one can write down the following standard model invariant 
1 Color sextet fields are preferred over color triplets, since the sextets do not mediate proton decay. 



4 



interaction Lagrangian: 2 

d = %x« + + + u c ad + ^sx 2 r + ^sxz 2 + h.c. (i) 

z J 2 J 2 J J 2 z 

If the scalar field S which has B = 2 is given a vacuum expectation value, cubic scalar field 
couplings of the type X 2 Y that break baryon number by two units will be induced. In turn 
it will lead to iV «-> N oscillation via the diagram of Fig. 1 with S r replaced by (S) [H] . 

We note that not all of the (X, Y, Z) fields are needed for B violation and N <-» N 
oscillation. (X, Y) or (X, Z) fields will do. We will focus more on these minimal versions in 
our computation of baryon asymmetry, while for generality we keep all three fields. 

To see the constraints on the parameters of the theory, we note that the present limits 
on r A r_jv — 10 8 sec. implies that the strength G N _fj of the AB = 2 transition is < 10~ 28 
GeV -5 . From Fig. 1, we conclude that 

r _ „ Ai (S) hjj u A 2 (S) hng^ 28 5 

For A1.2 ~ 1, hn ~ fix ~ gn ~ 10 -3 , we find (S) ~ Mxy.z — 1 TeV is allowed. In our 
discussion, we will stay close to this range of parameters and see how one can understand 
the baryon asymmetry of the universe. In fact, we will see that the masses of X, Y, Z cannot 
be much larger than a TeV for successful baryogenesis. Note that the couplings (/, g, h)ij to 
the second and third generation fermions could be larger. 

Other constraints can come from low energy observations such as bounds on flavor chang- 
ing hadronic processes such as K — K, D — D transition etc. This of course depends on any 
possible mixings between the right handed quark fields, on which we do not have any apriori 
information. If we make the simplest assumption dictated by the left-right symmetric theo- 
ries that the left and the right-handed mixings are equal, then the strongest constraints come 

2 An additional term g'^ Z (u^d^ — u^df) with g'^ = —g'^ is also allowed by the standard model gauge 
symmetry, but this term is forbidden when the model is embedded minimally into a left-right symmetric 
framework. We do not keep this term explicitly here, its inclusion is however straightforward. 



from K — K transition which imply that for h\\ ~ 10 3 , Mx > 1 TeV, which is consistent 
with our choice of parameters dictated by observability of iV <-> N transition. 

The model of Eq. (1) is embeddable into an SU{2) L x SU{2) R x S77(4) c framework 
where the quarks and leptons transform as ip : (2,1,4) © %p c : (1,2,4) representations 
and the Higgs fields X, Y, Z, S are part of the A c : (1, 3, 10) multiplet. In fact the S field 
corresponds to the A£ Cj , c component that acquires a vev and breaks B — L by two units. In 
this paper, we will not discuss the full set of constraints that arises in this embedding but 
rather simply work within the scalar field model described in Eq. (1). All our conclusions 
below apply to the SU{2) l x SU(2)r x SU (4) c model as well. While we take baryon number 
as part of the gauge symmetry, the mechanism of B asymmetry generation also works if B 
is a spontaneously broken global symmetry as in Ref. [T3| . 

3 Origin of matter 

Before proceeding to the discussion of how baryon asymmetry arises in this model, let us 
first consider the effect of the new interactions in Eq. (1) on any pre-existing baryon asym- 
metry. For this purpose, we assume the following mass hierarchy between the S field and 
the (X, Y, Z) fields: M s ~ 100 GeV < M x ,y,z ~ TeV. For T > M x ,y,z, the AB = 2 in- 
teraction rates scale like T and are in equilibrium at least down to T ~ Mx,y,z- They will 
therefore erase any pre-existing baryon asymmetry. They remain in equilibrium down to the 
temperature T* determined by the inequality: 

(2tt)9 M| 2 z ~ M Pl ■ U 
Here h and g refer to the largest of hij or gij are family indices). For h, g in the range 
of 0.1 - 1, this leads to T* ~ (0.6 - 0.2)M x ,z- 

The singlet field S will play a key role in the generation of baryon asymmetry. We assume 
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that (S) ~ Mx and M Sr ~ 10 2 GeV, where S r is the real part of the S field after its vev is 
subtracted. S r can then decay into final states with B = ±2, viz., S r — > 6q and SV — > 6g, 
inducing a net baryon asymmetry. 

On the way to calculating the baryon asymmetry, let us first discuss the out of equilibrium 
condition. As the temperature of the universe falls below the masses of the X, Y, Z particles, 
the annihilation processes XX — > d c d c (and analogous processes for Y and Z) remain in 
equilibrium. As a result, the number density of X, Y, Z particles gets depleted and only 
the S particle survives along with the usual standard model particles. The primary decay 
modes of S r are S r — > u c d c d c u c d c d c and S r — > u c d c d c u c d c d c . There could be other decay 
modes which depend on the details of the model. Those can be made negligible by choice 
of parameters which do not affect our discussions of N <-> N oscillation and baryogenesis. 
We will discuss them later in the paper. For T > M Sr) the decay rate of S r is given by the 
left-hand side of Eq. (3). This decay goes out of equilibrium around T* ~ 0AM x , or around 
500 GeV. Below this temperature the decay rate of S r falls very rapidly as the temperature 
cools. However as soon as T < M Sr , the decay rate becomes a constant while the expansion 
rate of the universe slows down. So at a temperature T d) S will start to decay where T d is 
given by 

T lSP\lh 2 g 4 M P ,Mf 1 1/2 
.(2tt) 9 1. QQgl /2 (6M x ) 12 _ 

This is obtained by equating the decay rate of S r to the expansion rate of the universe. In Eq. 
(4) the factor 18 is a color factor, h 2 = Tr(hJh), etc, while P is a phase space factor, which we 
have computed for the six body decay via Monte Carlo methods and found P ~ 2.05. The 
corresponding epoch must be above that of big bang nucleosynthesis. This puts a constraint 
on the parameters of the model. For instance, for M s ~ 200 GeV and M x ~ TeV, we get 
T d ~ 40 MeV (for g ~ h ~ 1). We will conduct the rest of the discussion with this set 
of parameters as a representative set. Note that M x cannot be much larger than about 1 
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Figure 1: Tree level diagrams contributing to S r decays into 6 anti-quarks. There are other 
diagrams where S r decays into 6 quarks, obtained from the above by reversing the arrows of 
the quark fields. 

TeV, otherwise T d will be below few MeV, affecting big bang nucleosynthesis significantly. 
Note also that the at least some of the couplings in h and g should be of order one. This 
would imply that the first family couplings should be of order (10~ 3 — f 0~ 4 ) from naturalness 
(tin ~ ^td^33 e ^ c )> making N <-> N oscillation accessible to next generation experiments. 

We now proceed to calculate the baryon asymmetry in this model. It is well known 
that baryon asymmetry can arise only via the interference of a tree diagram with a one 
loop diagram which has an absorptive part. The tree diagrams are clearly the one where 
S r — > 6q and S r — > 6q. There are however two classes of loop diagrams that can contribute 
to baryon asymmetry: one where the loop involves the same fields X, Y and Z as in Fig. 2 
(there is a second loop diagram of this type with (X, Z) fields inside the loop), and a second 
one involving W ± gauge boson exchange as shown in Fig. 3. In the (X, Z) model and in 
the (X,Y) model, only the latter contribution exists (the former trace being real). So we 
focus on that latter, which involves only standard model physics at this scale and has the 
advantage that it involves less number of arbitrary parameters. 
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Figure 2: Diamond loop diagram in the (X, Y, Z) model. 

We summarize the results of our calculations for the exchange diagrams. If one 
of the external up-type quarks is the top quark, the corresponding quark line receives a 
wave function correction via W gauge boson exchange. The baryon asymmetry from this 
diagram is found to be 



- ■ (5) 

Br 8 \^ ' mf J m t rn^{gg^) m 

where M u = diag(m u , m c , m t ), M& = diaglma, m s , m^) and V is the CKM matrix. Br 
stands for the branching ratio of S r into 6q + 6q. 

The vertex correction via the W boson exchange gives an asymmetry given by 

e£ ertex a 2 \uiYi[g T M U V gW* M d ] 



(6) 



Br 4 Tr(g*g) 

Here we have assumed that Ms r 3> m t . In the limit where rris r -C mw, we have the same 

asymmetry as in Eq. (6), but with a factor of (-1/4) multiplying it. Of course in this case, 

decays involving final state top quark are disallowed, which is to be implemented by removing 

the top quark contribution in the trace of Eq. (6). 3 

3 These loops do not conflict with the theorem of Ref. which states that no baryon asymmetry 




Figure 3: One loop vertex correction diagram for the 5-violating decay S r — > Qq. There are 
also wave function corrections involving the exchange of gauge boson. 

It is interesting to note that in this mechanism, there is a natural explanation of the 
observed baryon asymmetry r\ B ~ 10~ 10 . It follows from the light quark mass and mixing 
angle suppression. As an example, consider the following choice of parameters: m c (m c ) = 
1.27 GeV, m b (m b ) = 4.25 GeV, m t = 17 A GeV, V cb ~ 0.04, M s = 200 GeV and 
Is^l — \923\ ~ 1) w hh smaller values of g u . We find e B ~ 10~ 8 in this case from Eq. 
(5). The corresponding value from Eq. (6) is an order of magnitude larger, for the same 
input parameters. 

There is a further dilution of the baryon asymmetry arising from the fact that T d <C 
M$ r since the decay of S r also releases entropy into the universe. In this case the baryon 
asymmetry reads 



In order that this dilution effect is not excessive, there must be a lower limit on the ratio 
T d /Ms r . From our estimate above we require that T d /M Sr > 0.01. Since the decay rate of 
the S r boson depends inversely as a high power of M^y, higher X,Y bosons would imply 

that T s ~ H is satisfied at a lower temperature and hence give a lower T d /M Sr . In Fig. 4 

can be induced by dressing a AB = 1 vertex by baryon number conserving interactions. Since S r field has 
no definite baryon number, owing to {S) ^ 0, the theorem is not applicable in our case. 



(7) 
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Figure 4: The allowed range of Mx and Ms needed to generate the baryon asymmetry (along 
the black curve), decay temperature above 200 MeV (points below the dashed curve) and 
t~nn > 10 8 sec with A = 10 -4 (points above the red curve). 

we have plotted Mxy vs M$ r which gives the right amount of baryon asymmetry which is 
consistent with the demand that the decay of S r occurs before the QCD phase transition (i.e 
T d > 200 MeV). Using the effective coupling A = (Ai^/n) 1 / 4 ~ (X 2 hgf 1 ) 1 ^ to be of order 
10~ 4 implies that the 10 9 sec < t n _n < 10 11 sec for M$ r in range of ~ 100 — 300 GeV . 

We now comment briefly on some other aspects of the model: 

(i) If the (X, Y, Z) scalars are all present, the loop diagram of Fig. 2 will contribute to 
baryon asymmetry. Since there are two diagrams of the type shown in Fig. 2 the relevant 
trace has an imaginary piece. The asymmetry will have a suppression factor Mg r /M x , in 
addition to a loop factor and the Yukawa suppression. Although not very predictive, this 
model can yield adequate baryon asymmetry. 

(ii) The singlet field S can have a renormalizable interaction with the standard model 
Higgs doublet of the form XsS^SH^H. After the fields S and H acquire vacuum expectation 
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values, the Re S and the SM Higgs field can mix with each other opening up new decay 
channels for the S r field such as S r — > bb etc. We estimate this decay width to be 

3(3 2 m 2 b M s 



T(S r -> bb) 



where (3 is the S r — h mixing angle. There are two constraints on this decay mode in order for 
our scenario to work: first, this could contribute dominantly to the S r decay width thereby 
diluting the baryon asymmetry. Second, this model should be out of equilibrium at T = M$. 
If the model is non-supersymmetric, these two conditions are satisfied if As < 1CT 6 . This 
coupling is automatically forbidden if the model is supersymmetric. 

(iii) The present considerations could be easily extended to include supersymmetry as 
part of the new physics beyond the standard model. The SX 2 Y and SXZZ interactions in 
this case are nonrenormalizable ^3] • However, in this case we also expect mass terms in the 
superpotential of the form MsSS so that the effective four scalar interaction responsible for 
baryogenesis is in the same form as discussed above. 

(iv) Our theory is also testable in collider experiments such as the LHC since we have 
colored diquark scalar fields with masses in the TeV range. In a pp collision one could 
produce the X, Y, Z bosons either in pairs via the process qq — > XX or singly via the 
process q + g — > X + q. In the first case the signal would be a four jet final state whereas 
in the second case, it would be three jet final states. It would therefore be important to 
search for such final states at LHC. One distinguishing feature of these bosons is that they 
carry baryon number, which may be testable in the decays of these bosons into top quark 
and bottom quarks. 
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4 Conclusion 



In conclusion, we have presented a new mechanism for baryogenesis which does not rely on 
the electroweak sphalerons but rather directly produces matter-antimatter asymmetry using 
higher dimensional baryon violating couplings of a scalar field. The baryon asymmetry is 
produced at the weak scale. This mechanism can be tested by searches for baryon violating 
processes such as neutron antineutron oscillation, as well as by the discovery of colored 
scalars at the LHC 
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